Layer-resolved conductivities in multilayer graphenes 
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We study interlayer transport of multilayer graphenes in magnetic field with various stacking 
structures (AB, ABC, and A A types) by calculating the Hall and longitudinal conductivities as 
functions of Fermi energy. Their behavior depends strongly on the stacking structures and selection 
of the layers. The Hall conductivity between different layers is no longer quantized. Moreover, for AB 
stacking, the interlayer conductivity vanishes around zero energy with increasing layer separation, 
and shows negative values in particular cases. The fact that longitudinal interlayer conductivity 
suppressed by the magnetic field indicates that this system can be applied as a switching device. 
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— Introduction — Graphene has attracted increasing 
attention since its first isolation from graphite in 2004 [l[ . 
Graphene consists of two-dimensional hexagonal lattice 
of carbon atoms, whose quasiparticles are governed by a 
massless Dirac equation. A variety of unusual phenom- 
ena are observed in this system such as universal value 
of the minimum conductivity, anomalous quantum Hall 
effect, and so on [l-Q- I n addition, bilayer graphene has 
also been studied intensely which is characterized 

by intrinsic Landau level degeneracy at zero energy and 
a gate tunable band gap. 

Due to the qualitative differences between mono- and 
bilayer graphene, a lot of interest has been focused on 
multilayer graphenes to determine, how additional layers 
influence their physical response. One of the most in- 
triguing property of these systems is the variety of stack- 
ing structures such as AB (Bernal), ABC (rhombohedral) 
and A A (simple hexagonal) types. Graphene is usually 
produced by micro- mechanical cleavage of graphite, so 
that the stacking structure is usually of AB type, since 
the natural graphite falls into this category. However, 
production of graphene with other stacking types is also 
possible by recent technology such as epitaxial meth- 
ods 0, [loj]. In addition, AA stacking can be realized 



by folding of a graphene sheet [ll|. In terms of band 
structure, multilayer graphenes with more than 10 lay- 
ers are regarded as bulk graphite [l2| . so that few layer 
graphenes have been considered to be important systems, 
interpolating between graphene and graphite. So far, dia- 
magnetism transport properties 14-2l[, and energy 
spectra 22|, [23| have been studied for these systems. 

In this Letter, we investigate the transport properties 
between two different layers of multilayer graphenes in 
magnetic field as illustrated in Fig. HJ using the Kubo 
formula. Here, while electric current is run through a 
given layer, the resulting voltage drop or induced cur- 
rent is measured in another layer. In order to calcu- 
late these "layer-resolved conductivities", we establish 
the formalisms to obtain eigenvalues and eigenstates of 
multilayer systems, using the block diagonalization tech- 
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FIG. 1: Schematic illustration of the layer- resolved mag- 
neto transport between /c-th and Z-th layers in a multilayer 
graphene. 



nique. We find many interesting properties such as neg- 
ative response and switching effect by magnetic field, ab- 
sent in monolayer graphene. 

— Formalism — We consider three types of the stack- 
ing structures of multilayer graphenes, AB, ABC and A A 
types. Since a monolayer graphene consists of two sub- 
lattices labeled by A and B, carbon atoms in multilayer 
graphenes are specified by (A^B^) meaning the A, B sub- 
lattices in the i-th layer, respectively. Figure [2] shows the 
lattice structure of these systems with nearest neighbor 
intralayer (interlayer) coupling t (t±). 

Among these three types, we discuss AB stacked 
graphene in detail which is the most common structure 
of graphite. After taking the continuum limit of tight- 
binding model, in a basis with atomic components for 
N layers, |Ai), |Bi), • • • , |Ajv), |Bjv), the model Hamilto- 
nian in the vicinity of K point (per spin and per valley) 



H V 
W H V f 
V Ho V 



^ 

ooo'-. '■ 



(1) 



with 



Ho 



V7T- 
V7T+ 



V 





t± 



(2) 



where 7r± = 7r x ± with 7r = p + eA/c being the mo- 
mentum operator in a magnetic field V x A = (0, 0, B). 
v = (y/3/2)at/h is the Fermi velocity with a being the 
lattice constant. We have ignored long-range hopping 
terms except for t and t± for simplicity. Since the com- 
mutation relation between the momentum operators in 
Eq. ([2]) is [tt± , 7Tzp] = =f2eBh/c, there are correspondences 
with the creation and annihilation operators of the har- 
monic oscillator: tt± — >• y/2j(rf and 7r : 
eB ^ 0, where I = 



V^f a for 



In order to solve this model, we employ the matrix 
decompositions of the Hamiltonian. It is already known 
that the Hamiltonians of AB and AA-stacked A^-layer 
graphenes can be block diagonalized, considering Fourier 
modes of the wave function along the stacking direction 
[ll, EMEU- The same conclusion can also be obtained 
by factorization of determinant of the Hamiltonians [l4| . 
According to these, the effective Hamiltonian of AB- 
stacked Af-layer graphenes can be divided into isolated 
[N/2]g effective bilayer systems ([x]g is the integer part 
of x), and one monolayer system is added when N is odd. 
Similarly, the effective Hamiltonian of an AA-stacked Al- 
layer graphene consists of N isolated monolayer systems 
with different potential energies. 

Therefore, we can introduce a transformation matrix 
U for the AB stacked system which relates the wave 
function in the real space | Ai), |Bi), • • • , |Ajv), |Bjv) 
and that in the Fourier modes of stacking direction 

\4> ( t e r\ \^LT n) ), \4> ( NL°i d) ), \4> ( n-T\ l^ 3 Ven) )r ' ' , 
| j&n), \^b)- Here we have used the notation defined in 
Refs.ll3l andl23l. Then the Hamiltonian % is transformed 



into a block diagonalized form, 
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Here, U sah {m) is a Hamiltonian of a bilayer system with 
an effective hopping 
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with 
\ rn = 2 sin 



2{N + 1) J ' 



ra = AT- l,N-3,--- >0. 
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FIG. 2: Lattice structures of (a)AB, (b)ABC and (^AA- 
stacked trilayer graphenes, containing six sites in a unit cell. 
White and black circles denote carbon atoms which belong to 
A and B sublattices in each layer. 



Using the above block diagonalized Hamiltonian, we eas- 
ily obtain eigenvalues of AB-stacked graphenes based on 
the known results for monolayer and bilayer systems, re- 
placing the interlayer hopping t± by the effective one 

t±Xm, 



+ W(A m r) 4 + 2(2n + l)(A m r) 2 + 1 
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where r 
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t± and n denotes the Landau levels. The 



label \± = (si, S2) specifies the outer and the inner bands 
(51 = ±1), and positive and negative (52 = ±1) energies, 
respectively. The eigenstates of AB-stacked graphenes 
in basis of the real space |Ai), |Bi), • • • , |Ajv), \&n) are 
obtained from those of the subsystems H suh (m) and the 
transformation matrices U, written as 

I*™,m>= [#>-!> Pn» /„ 3 » tijn + l) ■■ 

(6) 

where 1 , denote coefficients of the wave function 
for the k-th layer, \n) is the number state of a, , and 
1 < /jl < 2N is band indices of the multilayer system. 
The conductivity is given by the Kubo formula as 



Re <Jij (O) = 



Im Uij(Q) 
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where with {i,j} G y} and is the Fourier 
transform of the current-current correlation function ob- 
tained after analytic continuation of the Matsubara form. 
The general expression of the conductivity of the multi- 
layer graphenes is obtained by extending the result for 
the bilayer graphene Q as 
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FIG. 3: Schematic illustration of intensity of the layer- resolved conductivity for AB-stacked multilayer graphene applying bias 
to the /c-th layer where k is (a) odd and (b) even layers counting along the direction of the magnetic field. The layer- resolved 
Hall ((Txy) and longitudinal (cr^) conductivities of 4- layer graphene in magnetic fields (B = 14T) for (c)AB, (d)ABC and (e)AA 
stacking. The bias is applied to the 1st layer. An inset in (c) shows a)^ x as functions of the magnetic field. 



with 

X(A,B;Sl) =Y,[(^n-A/h)- 1 (ioj n+m - B/h)-\ Vm ^ 
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where Cj n is Matsubara frequency of fermion including 
the chemical potential and effect of impurity scattering 
T as iu) n + [/i+isgn(co> n )r]/fL v m is the bosonic Matsubara 
frequency and cj n+m = uj n + v m - 

In the above formalism, the electric current operators 
are defined by Ji = — jj- with Ai the vector potential in 
the direction i = x,y. Now, we introduce current opera- 
tors for particular layers to calculate the conductivity be- 
tween two distinct layers as Jf = —j^k where k denotes 

layer number (k = 1, 2, • • • , TV), and is the vector po- 
tential acting solely in layer k in direction i. The current 
operator of the k-th layer has only two matrix elements, 
(Ji)2k-i,2k an d (Ji)2k,2k-i, an d all other elements are 
vanishing. The general expression for the layer-resolved 
conductivity between the k-th and l-th layers, afj is given 
by Eq. (|8|) with the following replacement, 
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In the present model which includes only the nearest in- 
ter layer couplings as inter layer matrix elements, the layer 
current operator, does not have any momentum de- 
pendence. However, this situation can change if we intro- 
duce other matrix elements, such as tilted inter layer hop- 
pings. In such cases, the layer current operators should 
be redefined appropriately so that they become Hermi- 
tian and satisfy the relation Ji = Ylk=i ^i- 

— Numerical results — Based on the above discussions, 
we calculate the layer resolved Hall conductivity a^y 
of AB-stacked multilayer graphenes, as functions of the 



Fermi energy, for strong magnetic field ftw c ^> T where 
uj c oc \fB is the cyclotron frequency. Our findings are 
summarized as: i) a^y 18 no longer quantized as inte- 
ger times e 2 /h. ii) around the zero energy fi ~ 0, a^y 
becomes almost zero when the source and the drain are 
separated by more than three layers \k — l\ > 3 and the 
biased layer k is odd, and two layers \k — l\ > 2 and k is 
even, iii) the layer resolved Hall response can be negative 
( a xy < 0) aroim cI the zero energy for \k — 1\ = 2 when the 
biased layer k is odd. Layers indices are growing in the 
direction of the magnetic field. These features are sum- 
marized in Fig. C3a),(b). Similar behavior can also be 
seen in the longitudinal conductivity a%x which behaves 
like derivatives of a^ l y by \i. 

In Fig. [3] we show numerical results of a% l y for TV = 4 as 
functions of the Fermi energy /i, where the bias is applied 
to the first layer, and B' = 2v 2 eB/c. We have assumed 
B = 14T, T = OK, T = 0.01*_l, t = 3.16eV, t± = 0.39eV, 
and taken DC limit ft — >> 0. We can see the above three 
features for AB staking in Fig. EJc). Further, is fi- 
nite without magnetic field, meaning that property ii) 
is broken when magnetic field is turned off, as shown in 
the inset of Fig. [31(c). The conductivity vanishes when 

b' ~ r. 

We also calculate layer resolved conductivities for other 
stacking structures. For ABC stacking, although the ma- 
trix decomposition technique is no longer available, we 
can diagonalize the Hamiltonian numerically to obtain 
the layer resolved conductivities. In this case, as shown 
in Fig. E](d), the electric currents are induced in every 
layer. Moreover, the conductivities for ABC-stacking do 
not exhibit negative values even in the vicinity of zero 
energy. 

For AA stacking, the matrix decomposition technique 
similar to AB stacking can also be used HI, [U[ . We cal- 
culate the conductivities by setting the value of interlayer 
hopping as t AA = t AB ' ABC /2. In this system, negative 
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conductivity also appears as shown in Fig. [3je). In this 
case, it is difficult to summarize the features of layer re- 
solved conductivities in terms of simple rules, because 
the effective Hamiltonian consists of N monolayer sys- 
tems with different Fermi energies, so that the number 
of Landau levels near the zero energy depends strongly 
on the strength of the magnetic field. In contrast to this, 
the Landau level structure is an intrinsic property for AB 
and ABC-stacked systems. 

— Analytical results for AB stacking — In order to un- 
derstand the above results for AB stacking in more detail, 
we calculate the first quantum Hall step of a^ l y analyt- 
ically for T = 0. We consider the following four cases: 
(a/b) (fc, Z) = (odd, odd) for N =even/odd, respectively, 
(c) (fc,/) = (odd, even), (d) (fc,/) = (even, even). Since 
layer indexes grow in the direction of the field, such classi- 
fication is unique. We obtain the following results, defin- 
sgn(e£)^<7^, 
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where fc' = [^r"]^ and Z' = From these results, 

it is apparent that interlayer conductivity is finite only 
within nearest or second nearest layers (property ii)). For 
the simplest example, the values of a%y for three layer 
system obtained from the above formalism are cr^y = 
(2 + 3r 2 + r 4 )/(4 + 8r 2 ), = r 2 /(4 + 8r 2 ), v% = 



a 1 ' 2 

xy 
4 / 



(2 + 2r 2 )/(4+8r 2 ), = -r 4 /(4+8r 2 ), with a% 3 
and a^ 3 = cr^ 1 . In this case, <j^ 3 becomes negative. By 
summing up all contributions, we obtain l a^. y l =3/2 
which consists of contributions from effective monolayer 
(1/2) and bilayer (1). For other N cases, these analytical 
results coincide with the general results summarized in 
Fig.EJa) and (b). 

— Conclusion and discussion — To conclude, we have 
studied the interlayer electronic transport properties of 
multilayer graphenes in a magnetic field for variety of 
stacking orders. The behaviour of the layer-resolved con- 
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ductivity depends strongly on the stacking structure. For 
AB stacking, various interesting properties, such as neg- 
ative response, and suppression of the Hall conductivity, 
are identified. The breakdown of the quantization of the 
Hall conductivity indicates that interlayer conductivity 
has different features from those of the total response, 
and is not protected by topology. 

Finally, we discuss the possibility of the experimental 
observation of the layer-resolved conductivities. Though 
it would be rather challenging to connect leads to partic- 
ular layers, it is certainly easier to measure the layer- 
resolved conductivity between the top and the bot- 
tom layers. The fact that the longitudinal conductivity 
greatly changes according to the magnetic field (for ex- 
ample shown in the inset of Fig. 02(c)) means that 
multilayer systems may be applied as switching device. 
We think that our work provides a comprehensive under- 
standing of transport properties of multilayer graphene. 
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